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We present a theoretical and experimental study of the damping process of the atomic velocity 
in Sisyphus cooling. The relaxation rates of the atomic kinetic temperature are determined for a 3D 
lin_Llin optical lattice. We find that the damping rates of the atomic temperature depend linearly 
on the optical pumping rate, for a given depth of the potential wells. This is at variance with the 
behavior of the friction coefficient as calculated from the spatial diffusion coefficients within a model 
of Brownian motion. The origin of this different behavior is identified by distinguishing the role of 
the trapped and traveling atoms. 



I. INTRODUCTION 

The last decades have witnessed an impressive progress in laser cooling techniques, and nowdays it is possible to 
prepare and manipulate very cold and dense atomic samples [1,2]. Sisyphus cooling [3] represents a milestone in the 
history of laser cooling, particularly because this new cooling process demonstrated that the Doppler limit could be 
broken [4]. Furthermore, it is an important peculiarity of the Sisyphus process that the cooling mechanism not only 
leads to very low temperatures, but also to the strong localization of the atoms near the minima of the optical potential 
associated to the dipole force. Using appropriate laser configurations, it is then possible to design multidimensional 
modulations of the optical potential such that the atoms result spatially ordered in periodic arrays whose periodicity 
is of the order of the laser wavelength, the so-called optical lattices [5]. The possibility to prepare atoms in spatially 
ordered structures turned out to be essential for several studies in nonlinear dynamics [6-8]. 

Since the early works [3], Sisyphus cooling has been extensively studied, both experimentally and theoretically, in 
monodimensional as well as multidimensional optical lattices. The dependence of the steady-state kinetic temperature 
on the lattice parameters has been carefully examined [9-11] . It is by now a well established result that the temperature 
scales as the depth of the optical potential [9] . More recently, it has been shown that in a 3D lin_Llin configuration (see 
below) the temperature is on one hand independent of the lattice constant [12,13] and on the other hand anisotropic 
[13,14]. Much less attention was payed to the study of the time scale of the damping of the atomic velocity in an 
optical lattice. So far, only indirect measurements of the cooling rate of atoms have been presented [15]. In that work, 
the cooling rates were inferred from the atomic localization rates measured by using Bragg scattering as a sensitive 
probe. 

In this paper we report a theoretical and experimental study of the damping process of the atomic velocity in 
Sisyphus cooling. More precisely, we measure, both in the experiments and in the numerical simulations, the width of 
the atomic velocity distribution (which is proportional to the square root of the kinetic energy of the atomic sample) 
for different lattice durations within the thermalization phase. It has thus been possible to monitor the damping 
of the kinetic energy in a direct way. We show that the kinetic energy decreases exponentially during this phase 
and we determine the corresponding damping rate Tt as a function of the lattice parameters (detuning from atomic 
resonance and intensity of the laser beams). We find that F^ is proportional to the optical pumping rate and is 
strongly anisotropic. We also compare the cooling rate F^ to the friction coefficient a as deduced from the spatial 
diffusion coefficients by using the Einstein relation for Brownian motion. We find that these two quantities show 
completely different dependencies on the lattice parameters. The origin of this different behavior is identified by 
distinguishing the role of the trapped and traveling atoms. 

The present investigation is motivated by a parallel study of probe-transmission spectroscopy of atoms cooled in 
a dissipative optical lattice [16,17]. The properties of the optical lattice are studied by monitoring the transmission 
of a weak probe field through the atomic sample as a function of the relative detuning between the probe and the 
lattice fields. The knowledge of the relaxation rates of the different atomic observables (atomic density, atomic velocity, 
magnetization, etc) is necessary to understand the physical origin of the Rayleigh resonance that appears at the center 
of the probe transmission spectra. It is well established that the Rayleigh resonance originates from the scattering 
of photons from the lattice beams on the grating of an atomic observable induced by the pattern of interference 
between the probe and the lattice beams, and that the width of the Rayleigh line is given by the relaxation rate of 
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this observable [18] . However which atomic observable determines the Rayleigh line in an optical lattice is still an 
open question [17]. Hence, the knowledge of the relaxation rates of the different atomic observables should allow, 
through comparison with the width of the Rayleigh line, the identification of the observable leading to the Rayleigh 
resonance. 



II. EXPERIMENTAL SET-UP 



In our experiment *^Rb atoms are cooled and trapped in a 3D dissipative optical lattice red detuned with respect 
to the Fg = 3 ^ Ff. = 4: D2 line transition. The procedure to load the atoms in the optical lattice is the standard one 
used in previous experiments [12]: the rubidium atoms are first cooled and trapped in a magneto-optical trap (MOT) 
for about 1 s; the temperature of the atoms in the MOT is of the order of 200 /iK. Then the MOT magnetic field 
and laser beams are turned off and the lattice beams are turned on by using acousto-optical modulators (AOM) with 
a fall-rise time shorter than 1 fis. 

The three-dimensional periodic optical lattice is generated by the interference of four linearly polarized laser beams, 
arranged in the 3D lin_Llin configuration [5] (Fig. 1): two j/-polarized beams propagate in the (xOz) plane with a 
relative angle 26, while two x-polarized beams propagate in the (yOz) plane and form also an angle 20. This 
arrangement results in a periodic modulation of the light polarization and light intensity, which produces a periodic 
modulation of the light shifts of the different atomic ground-state sublevels (optical potentials). The optical potential 
has minima located on a orthorombic lattice and associated with pure circular (alternatively cr+ and a~) polarization. 
The lattice constants, i.e. the distance (along a major axis) between two sites of equal circular polarization are 
^x,y = X/sinO and = A/(2cos(?), where A is the laser field wavelength. 

This configuration is well known to be suited for Sisyphus cooling [5] for which the steady-state kinetic energy of 
the atomic sample results from the competition of cooling and heating processes. The optical pumping between the 
different sublevels of the atomic ground state combined with the spatial modulation of the light shifts leads to the 
cooling of atoms, while the heating is esscntialy due to the statistical fluctuations of the dipole force [3]. For deep 
enough potential modulations, Sisyphus cooling leads to atomic kinetic energies lower than the potential depth so 
that the atoms get trapped at the minima of the optical potential wells, thus producing a periodic array of atoms. 

The temperature of the atoms is measured with the method of recoil-induced resonances [19-22]. The atoms are 
released from the lattice by turning off the lattice beams and two new laser beams (the pump (c) and the probe (p) 
beams, see Fig. 1) are introduced. Beams c and p are detuned by 49 MHz from atomic resonance and their intensities 
are of the order of 100 /uW/cm^. They cross the atomic sample in the (xOz) plane, and they are symmetrically 
displaced with respect to the z— axis forming an angle of 2(p = 25°. The diamater of the pump beam is about 1 cm, 
larger than the atomic cloud, while the probe beam is very thin (less than 500 /im in diameter). The frequency of the 
probe beam is scanned around the pump frequency in 1 ms. We notice that for these parameters of the pump and 
probe beams, the heating due to their absorption is negligible. The probe transmission is monitored as a function 
of the detuning S between the pump and probe fields. For equal pump and probe polarizations, the ground states 
\a,p) and \a,p') with the same internal quantum number a but different atomic momenta are coupled by Raman 
transitions. The difference in population of these states determines the absorption or gain of the probe, and the 
probe transmission spectrum is thus proportional to the derivative of the atomic velocity distribution [20]. These 
distributions were found to fit well to Gaussian functions so that it is straightforward to calculate the mean kinetic 
energy Ek in the j- direction and to derive an atomic kinetic temperature 

Ek^ = ^ (1) 

from the width of the resonance in the probe transmission spectrum. The geometry of our pump and probe fields 
corresponds to measurement of the temperature in the x— direction. 



III. EXPERIMENTAL RESULTS 



In the experiment, the duration of the optical lattice is carefully controlled and can be varied in steps of 10/xs. We 
measure the kinetic temperature of the atomic cloud as a function of the duration t of the optical lattice for different 
lattice parameters (intensity per lattice beam and lattice detuning A from atomic resonance). We plot in Fig. 2 a 
typical variation of the measured kinetic temperature Tx as a function of the time t together with a fit of the data to 
a function of the form 

Tx it) = Tf^ + m - Tf^ ) exp(-rT. • t) . (2) 
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The kinetic temperature decreases exponentially with the lattice duration t towards a steady-state temperature, and 
the fit to Eq. (2) provides the steady-state temperature Tj^ as well as the damping rate of the atomic temperature. 
From experimental data as those of Fig. 2, we determined Tt.^ and T/^ in a broad range of lattice parameters (intensity 
per lattice beam II and lattice detuning A). We found that the temperature shows a linear increase with II and 
decreases against A, in agreement with earlier studies. A more complete account of the dependence of the steady-state 
temperature on the lattice parameters has already been reported in Ref. [12] and will not be repeated here. The novel 
result of the present work consists instead in the determination of the damping rate of the temperature. 

Experimental results for Tt^ are reported in Fig. 3 as a function of the lattice beam intensity at fixed lattice 
detuning (top) and as a fimction of the optical pumping rate Fq oc /^/A^ at fixed light shift per beam Aq oc /l/A 
(bottom). We find that F^^ ranges from a few kHz to a few tens of kHz for typical intensity and detuning of the lattice 
beams. More precisely, from the experimental data, we identify the following general behavior of the damping rate as 
a function of the lattice parameters: Ft-^ is an increasing fimction of the lattice beams intensity 7^, the dependence 
being linear within the experimental error, and a decreasing function of the lattice detuning |A|. When plotted against 
the optical pumping rate, Ft^ is found to be proportional to Fq, and almost independent of the light shift per beam. 
The observed behavior is in good agreement with the results of Raithel et al [15], and will be interpreted with the 
help of Monte-Carlo simulations. 

IV. THEORETICAL ANALYSIS 

We analyze theoretically the damping process of the atomic velocity with the help of semi-classical Monte-Carlo 
simulations. Taking advantage of the symmetry between the x— and y— directions (see Fig. 1) we restricted the atomic 
dynamics in the (xOz) plane. The calculations are made for a Jg = 1/2 ^ Je = 3/2 atomic transition, as customary 
done in numerical analysis of Sisyphus cooling. The numerical calculations in this simplified configuration are expected 
to give the same dependencies of the dynamical quantities on the lattice parameters as full 3D calculations, so they 
are useful to study the main fcauturcs of the cooling dynamics. However, a direct quantitative comparison to the 
experimental results is not necessarily meaningful because the restriction of the dynamics to two dimensions may 
introduce scaling factors [13]. 

A. Cooling rates 

We first determine the mean kinetic energies along the x— and z— directions as functions of time in the thermal- 
ization phase. The mean kinetic energy is obtained by averaging over about 5000 trajectories. For both directions, 
the temperature shows an exponential decrease with the lattice duration. By fitting the temperature and with 
a function of the form (2), we determine the corresponding relaxation rates Fy^ in the x— and z— directions, with 
results as in Fig. 4. 

We notice that the damping rates for the two orthogonal directions are very different, with Ft^ typically about 
one order of magnitude smaller than F^^ . This property certainly relates to the strong anisotropies of the optical 
potential and of the optical pumping rates in the 3D lin_Llin lattice. Previous work already pointed out anisotropic 
temperatures [13,14] and anisotropic spatial diffusion coefficients [12,13] for a 3D-lin_Llin lattice. 

The graphs of the left-hand column of Fig. 4 report the results for the damping rates Tt^ as functions of the intensity 
per lattice beam II for different values of the lattice detuning A and confirm the general behavior observed in the 
experiment: Tt^ is an increasing function of the lattice intensity, and a decreasing function of the lattice detuning. 

On the right-hand column of Fig. 4 the numerical results for the damping rates are reported as functions of the 
optical pumping rate per lattice beam Fq for fixed values of the light shift per beam Aq, i.e. for fixed values of the 
depth of the potential wells. This allows a more transparent interpretation of our results. Indeed, the damping rates 
are foimd to be proportional to the optical pumping rate in a large range of light shifts per lattice beam within a 
numerical error of about 10 %. This is in qualitative agreement with our experimental observations, and also with 
the indirect measurements of the cooling rates of Raithel et al [15]. 

B. Priction coefficients 

It is interesting to compare the results obtained for the damping rate of the kinetic energy with the well known 
friction coefficients for Sisyphus cooling [3]. Such a coefficient a is introduced in the low- velocity limit where the 
mean cooling force has the form F = —av. We determined the friction coefficients and for the two orthogonal 
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directions. A possible way to calculate the friction coefRcients is to introduce an arbitrary constant force Fj along 
the direction j. The mean friction force compensates this force, with the atomic cloud set into a motion along j at 
a constant mean velocity (vj) = Fj/uj. The friction coefficient aj is then easily derived by calculating the velocity 
{vj) for various Fj. A different approach to the calculation of the friction coefficients is based on the assumption that 
the atomic dynamics can be described as a Brownian motion. Under this assumption, the friction coefficients can 
be derived from the steady-state kinetic temperature Tj and the spatial diffusion coefficient Ds^ in the j-direction 
by using the Einstein relation for Brownian motion aj = kBTj/Dgj. In previous work [13], these two approaches 
have been compared and shown to lead to the same values for the friction coefficients. Also for the present choice of 
interaction parameters, we found that the two methods lead to values of aj in excellent agreement [23]. These values 
are shown in Fig. 5 as functions of Fq for different values of the light shift per beam Aq. We notice that also for the 
friction coefficients there is a strong anisotropy, with a^, about 5 times smaller than a^. The results for the friction 
coefficients have to be compared with the values of Ft^ and Tt^ shown in Fig. 4. The values aj/M are about one 
order of magnitude larger than those for ■ Furthermore, the two quantities do not show the same dependence on 
the lattice parameters. In fact, aj weakly depends on Fq but strongly depends on Ag in the range of parameters 
considered here. This behaviour is completely different from the one for Ft^ in the same range of lattice parameters. 

C. Discussion 

The very different dependence on the lattice parameters shown by the damping rates of the kinetic energy and 
the friction coefficients is rather surprising. In fact, the description of the atomic dynamics as a Brownian motion 
suggests that the friction coefficient in a given direction, rescaled by the atomic mass M, is simply the half of the 
damping rate of the kinetic energy in the same direction. The different behavior shown by these two quantities can 
be explained by distinguishing the role of the trapped and traveling atoms. In fact both in the thermalization phase 
and at thermal equilibrium there are high energy atoms (h-e) that travel over the lattice sites and low energy atoms 
(1-c) that are trapped in the potential wells. As discussed in detail in Ref. [13], the high energy atoms experience 
Sisyphus cooling, undergo a Brownian motion, and their dynamics can be characterized by a friction coefficient aj. 
Since these atoms undergo a Brownian motion one may expect 

On the other hand, low-energy atoms are trapped in the potential wells and do not undergo a Brownian motion. 

Therefore they show a very large (ideally divergent) friction coefficient. In fact, the application of a weak force Fj 
does not set them into a uniform motion, but only changes their center of oscillation in the wells. Therefore for these 
atoms the friction coefficient aj = Fj/{vj) is divergent. Furthermore, trapped atoms do not loose kinetic energy by 
Sisyphus cooling and therefore they are characterized by a cooling rate F^ ideally equal to zero. It follows that for 
the low-energy trapped atoms 




It is then clear that the presence of trapped atoms invahdates the simple relation aj/M = Ftj/2 and lead to a coohng 
rate smaller than 2aj/M, as observed in our numerical simulations. 

V. CONCLUSIONS 

In summary, we have determined experimentally the damping rates of the atomic kinetic energy (kinetic temper- 
ature) in Sisyphus cooling, and have examined their general dependencies on the lattice parameters. We have also 
made numerical semi-classical Monte-Carlo simulations. We have shown that the damping rates of the atomic kinetic 
temperature depend linearly on the optical pumping rate and are strongly anisotropic. We have also determined the 
friction coefficients from the atomic spatial diffusion cocfficincts and we have shown that they weakly depend on the 
optical pumping rate but strongly depend on the potential depth. The very different behavior of the damping rate 
of the kinetic temperature F^^ and of the friction coefficient Uj has been explained by distinguishing the role of the 
trapped and traveling atoms. 

The presented results for the damping rates of the atomic kinetic temperature in Sisyphus cooling are useful for 
the analysis of the Rayleigh resonance detected in probe-transmission spectroscopy of atoms cooled in dissipative 



4 



optical lattices [16,17]. Indeed, the knowledge of the damping rates of the different atomic observables excited in the 
optical process is necessary to determine the physical origin of the resonance. Preliminary measurements for the probe 
transmission show that the width of the Rayleigh resonance is of the same order of magnitude of the value for 
presented in this work. This suggests that the atomic velocity is one of the observables which, excited in the optical 
process determined by the probe and lattice fields, contribute to the Rayleigh resonance. Further work is however 
necessary for a definitive conclusion. 
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FIG. 1. Laser fields configuration for the 3D-lin_Llin optical lattice. The beams 1—4 generate the 3D periodic potential. For 
the measurement presented in this work, the angles between the lattice beams are kept fixed at ^ = 30°. Two additional laser 
beams (c and p), are introduced for the temperature measurement. They are both linearly polarized along the y-axis. The angle 
between beams c and p is 2<p = 25° . 
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FIG. 2. Experimental results for the temperature Tx in the x-direction of the atoms cooled in the optical lattice as a function 
of the duration t of the lattice. The lattice detuning is A = 48 MHz, the intensity per lattice beam II = 4.5 mW/cm^. The 
solid line is the fit to the exponentially decreasing function Eq. (2). 
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FIG. 3. Experimental results for the damping rate Ft^ of the atomic temperature in the x— direction. In the upper graph 
Ftj. is plotted as a function of the intensity per lattice beam II, for different values of the lattice detuning A. In the lower 
graph Ftj. is plotted as a function of the optical pumping rate Fq, for different values of the light shift per beam Aq. Here u>r 
is the atomic recoil frequency, and 7o the saturation intensity. 
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FIG. 4. Numerical results for the damping rates of the atomic temperature in the x— and z— directions. In the left hand 
column the damping rates are plotted as functions of the intensity per lattice beam II , for different values of the lattice detuning 
A. In the right hand column, the same data arc reported as functions of the optical pumping rate Fq, for different values of 
the light shift per beam Ag. Here F is the width of the excited state. 
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FIG. 5. Numerical results for the friction coefficients in the x— and z— directions as functions of the optical pumping rate 
To for various values of the light shift per beam Aq. 



